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Abstract. In this talk I will emphasize the role of the Truncated Overlap 
Fermions in showing the equivalence between the Domain Wall and Overlap 
Fermions up to an irrelevant factor in the fermionic integration measure. 
I will also show how Domain Wall type fermions with a finite number 
of flavors can be used to accelerate propagator calculations of their light 
partner in the infinite flavor limit. 



1. Introduction 

It required some time until Domain Wall [1, 2] and Overlap [3] formulations 
of chiral lattice fermions gained the adequate momentum [6]. A remnant 
chiral symmetry on the lattice, called the Ginsparg- Wilson relation [9], that 
was more recently noticed [4, 5], was shown to be the building block of a 
chiral gauge theory which exists on the lattice [7]. 

The basic idea of Domain Wall Fermions is an expanded flavor space 
which may be seen as an extra dimension with left and right handed 
fermions defined in the two opposite boundaries or walls, as it is sketched 
schematically below. 

Let N be the size of the extra dimension, Dw the Wilson-Dirac oper- 
ator, and m the bare fermion mass. Then, the theory with Domain Wall 
Fermions is defined by the action [1, 2]: 



N 

S DW := ^a 5 M DW ^ = ^2^i[(a B DW - 1)^ + P+^i+i + P-i/>i-i] (1) 
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with boundary conditions given by 

P+{ipN+i = () 

P-{ijj + mi/j N ) = 1 j 

where M. is the five-dimensional fermion matrix of the regularized theory 
and 1)11 = M — Z?vk with M G (0,2) being a mass parameter and 05 the 
lattice spacing in the 5 th direction. 




The theory with Truncated Overlap Fermions is defined by [8, 15]: 

N 

S DW := *a 5 M T ov* = Y,^ a 5 Dll - 1 )^iHa5D ll +l)P + iJi +1 +(a 5 D^+l)P^ i _ 1 ] 

i=i 

(3) 

In both cases the lattice spacing a of the four dimensional theory is set to 
one. 

Truncated Overlap Fermions can be formally constructed from Domain 
Wall Fermions by substituting 



P+Y> m -> (a 5 DW + l)P+A+i 
P-A-i (o 5 l>ll + 

while the boundary conditions remain the same as before. 

2. Continuum limit in the 5 th dimension 

Let me first write down the operator kernels of both theories: 

M DW = D\\ + ^(e a57595 - 1) 

Mtov = D\\(e as ^ d5 + 1) + ^{e abl6&5 - 1) 

This form can be easily checked by using the identity: 

^(*5 + a 5 ) = e a ^(t 5 ) 
with ipi = ip(ts = 10,5) 



(4) 



(5) 



(6) 
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Taking the limit as — > 0, I get: 

M D w = DW+ l5 d 5 + ^dl 

Mtov = DW{2 + a 575 c>5) + 75 5 5 + f 1 j 

Hence, in a continuous flavor space, both theories are unique up to 
an asymmetric factor remaining from the Truncated Overlap. Therefore, 
one may conclude that Domain Wall and Truncated Overlap Fermions are 
discretizations of the same Domain Wall fermion theory in the continuous 
flavor space [0, T5]: 
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defined by the following action: 

S = ^(D + l5 d 5 -M)W (8) 

with the following boundary conditions: 

P + [tf(-,T 5 )+m*(-,0)]=0 , , 

P_[tf(-,0)+mtf(-,r 5 )]=0 {) 



3. Ginsparg- Wilson relation 

A remnant chiral symmetry on the lattice may be possible if one allows a 
local symmetry breaking for propagating states. This statement is encoded 
in the Ginsparg- Wilson relation [9]: 

75D- 1 + LTSs = 2a 75 P, (10) 

where D is a local Dirac operator and R is also a local operator trivial in 
Dirac space. 

An explicit solution of this relation is given by the Overlap Dirac oper- 
ator [10]. In fact one can show that a Dirac operator obeying the Ginsparg- 
Wilson symmetry can be derived from the Domain Wall [11, 8], and Trun- 
cated Overlap Fermions [8] in the infinite flavor limit. 

The situation is unclear when the number of flavors is finite [8]. I present 
here some preliminary tests of the Ginsparg- Wilson relation on a small 
number of configurations on a 4 4 lattice at j3 = 6. 
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In Figs. 1-2 the locality of the Dirac operator is observed for N = 4 and 
N = 32 number of flavors. The behavior of R is tested in Figs. 3-4. These 
suggest that R tends towards a Kronecker-Delta function as the number 
of flavors grows and the convergence is faster for Domain Wall Fermions. 
More data are needed to verify this evidence. 



4. Infinitely separated walls 

The results of the previous section, although preliminary, are enough to 
conclude that the infinite limit in the fifth dimension is needed. This may 
be unrealistic for practical computations, if one would keep working with 
the whole 5— dimensional theory. 

A simple solution is to work in the four dimensional framework of the 
Overlap Dirac operator [10]: 

„ 1 + m 1 — m i \ , , 

D = — — 7ssgn(#) (11) 

where the Hamiltonian H = 75-D" corresponds to the "evolution" in the 
fifth dimension of Truncated Overlap Fermions with a transfer matrix given 
by [8]: 

Ttov = (12) 

For Domain Wall Fermions it is not straightforward to construct "easy 
to use" Hamiltonians, since the transfer matrix is given by [8]: 

T DW= 1 + ] Ip 0--HP+) (13) 

where numerator and denominator do not commute. 

In analogy to Truncated Overlap Fermions, I define a Hamiltonian H 
for Domain Wall Fermions, such that the transfer matrices of both theories 
coincide: 

from which I can write down the solution: 

D" 1 
"-»23Dir-*2TDir < 15 > 

where 05 = 1 is assumed. 

This looks merely a trick, but in fact it is obvious by the definition that 
H derives from the transfer matrix of the Domain Wall Fermions. Therefore, 
I arrive to the conclusion that 
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The light fermion operator in the infinite flavor limit of Domain Wall 
Fermions is given by the Overlap Dirac operator with Hamiltonian TL. 
Some remarks are in order here: 

a) The form of TL suggests that both theories are identical in the limit 
a — > 0. In this case TL « H. 

b) For finite a any theory with Wilson fermions can be equivalently 
defined to a theory with a Dirac operator: 

Pll _ D w 

2^m i+d w (16) 

for M = 1 and m = 
up to the determinant factor det(l+LV). This is easily seen by the identity: 

/ _ e -(x-^)(x-y>)-V>iVV> = det (! + D w )e~ XT ^ x (17) 
Jipip 

i.e. the new operator is the Schur complement of the new "effective" theory 
with free fermions x> X- Therefore up to an irrelevant determinant factor, 
both theories are equivalent for finite a. 
Computational remarks on TL. 

It is important to know the computation overhead of TL if one would 
like to work with Domain Wall Fermions in the infinite flavor limit, i.e. in 
the Overlap framework. 

Practical methods to compute the Overlap operator use the application 
of TL or TL 2 to a vector [12, 13, 14]. 

It is obvious that the computation of TL is more complex than that of 
H, although the inversion of 2 — is well conditioned and can be done 
fast. 

On the other hand, TL is conditioned better than H. To illustrate this, 
I have computed the spectrum of Z)"/(2 — I)") for free fermions on a 16 4 
lattice and also for a fixed background at j3 = 6 on a 4 4 lattice. The spectra 
are shown in Figs. 5-7. 

5. Inversion of the Overlap Dirac operator 

It has been pointed out that Truncated Overlap Fermions can be used to 
compute efficiently the inverse of the Overlap Dirac operator [8, 15]. From 
the discussion above, it can be concluded that Domain Wall Fermions can 
also be used effectively to compute the propagation of the light fermion in 
the infinite flavor limit. 

The basic idea is a multigrid algorithm, which is illustrated below. 



N = oo 

I 1 

A T 5 



| 1 1 1 1 1 1 1 1 

1 N 



I would like to solve the linear system: 

Dz = b (18) 

where D is the chiral Dirac operator and zq is a first guess. The algorithm 
may be described as a three step iteration scheme [8, 15]: 

I. Compute the quark propagator in the Domain Wall framework, i.e. 
finite N, which can be interpreted as a coarse lattice propagator: 

MPx 5 = MiPh 

6 5 = (6,0,...,0) T (19) 
X5 = (x,y,...,z) T 

where X5 and 65 are block- vectors with TV blocks, Mi is the same matrix 
M but with bare quark mass m = 1, and P is the following permutation 
operator: 

{Px 5 )i = P+(x 5 )i + P-(x 5 ) i+ i, i = l,...,N -1 , , 

(Px 5 ) N = P + (x 5 ) N + P_(x 5 ) 1 {ZU) 

II. Compute the residual error in the Overlap framework: 

Z = Zq + X 

r = b — Dz 



(21) 



III. Construct the new residual error of the five dimensional theory and 
define the new approximate solution: 

6 5 -(r,0,...,0f 

zo <- z 

and go to step I. or otherwise stop. 

The scheme is tested on 30 small 4 4 lattices at (3 = 6 for the Overlap 
Fermions. The results are shown in Fig. 8, where the multigrid pattern of 
the residual norm is clear. For comparison, in Fig. 8 are shown the results 
of directly applying the Conjugate Residuals (CR) algorithm. The gain is 
about a factor 10 in this case. More results are needed on larger lattices. 

Note that CR is the best, i.e. the optimal algorithm for the Overlap 
operator, which is a normal operator [16]. To invert the "big" matrix in 
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step I., I have used the BiCGstab2 algorithm [17] which is almost optimal 
in most of the cases for the non-normal matrices as its is the matrix M 
[16]. 



6. Conclusions 



I have shown the equivalence between Domain Wall and Overlap Fermions 
up to an irrelevant factor in the fermionic integration measure. 

Domain Wall and Truncated Overlap Fermions can be used to accelerate 
the computation of wall propagators in the infinite flavor limit. 
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Truncated Overlap Fermions 



o N = 4 
* N = 32 



Figure 1. Norm of the D kernel in spin and color space with the distance r from the 
origin for N = 4 (circles) and N — 32 (stars) . 



Domain Wall Fermions 




Figure 2. Norm of the D kernel in spin and color space with the distance r from the 
origin for N = 4 (circles) and N = 32 (stars) . 
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Truncated Overlap Fermions 
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Figure 3. Norm of the R kernel in spin and color space with the distance r from the 
origin for N = 4 (circles), N — 32 (stars) and iV = 64 (crosses). 



Domain Wall Fermions 
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Figure 4- Norm of the R kernel in spin and color space with the distance r from the 
origin for N = 4 (circles), iV = 32 (stars) and iV = 64 (crosses). 
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Free Fermions M = 1 16' 
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Figure 5. The spectrum of the D"/(2 — D") matrix in the complex plane. 
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Figure 6. The spectrum of the Z?"/(2 — D") matrix in the complex plane. 
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Figure 7. The spectrum of the Z)"/(2 — £>") matrix in the complex plane. 
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D w multiplications x1Q 4 

Figure 8. Norm of the residual error vs. the number of Dw multiplications on 30 
configurations. Circles stand for the straightforward inversion with CR and stars for the 
multigrid algorithm. 



